Intramolecular charge ordering in the multi molecular orbital system (TTM-TTP)l3 
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Starting from the structure of the (TTM-TTP)l3 molecular-based material, we examine the characteristics 
of frontier molecular orbitals using ab initio (CASSCF/CASPT2) configurations interaction calculations. It is 
shown that the singly-occupied and second-highest-occupied molecular orbitals are close to each other, i.e., this 
compound should be regarded as a two-orbital system. By dividing virtually the [TTM-TTP] molecule into 
three fragments, an effective model is constructed to rationalize the origin of this picture. In order to investigate 
the low-temperature symmetry breaking experimentally observed in the crystal, the electronic distribution in a 
pair of [TTM-TTP] molecules is analyzed from CASPT2 calculations. Our inspection supports and explains the 
speculated intramolecular charge ordering which is likely to give rise to low-energy magnetic properties. 



I. INTRODUCTION 

The characterization of molecule-based electronic conduc- 
tors is one of the central issue in the research of molec- 
ular crystal systems. Typical materials are TTF-TCNQ 
and (TMTSF)2PF6, where TTF, TCNQ, and TMTSF stand 
for tetrathiafulvalene, tetracyanoquinodimethane, and tetra- 
methyltetraselenafulvalene, respectively. In (TMTSF)2PF6, 
the first superconductivity behavior in molecular solids has 
been reported.^ Recent hot topics in the research of molecu- 
lar crystals is the realization of charge ordering phenomena.^ 
Since the pioneer work of Su et al? in poly acetylene, the 
charge trapping phenomenon has been much studied in Peierls 
transition issues.^ In that sense, quasi-one-dimensional chains 
have received much attention from experimental and theoret- 
ical points of view.^~^ Theoretical approaches that focus on a 
single highest-occupied-molecular orbital (HOMO) or lowest- 
unoccupied-molecular orbital (LUMO) have been successful 
in describing fascinating electronic ordered phase. Such a 
treatment can be justified since in these conventional systems, 
the HOMO or LUMO levels are well- separated from the rest 
of the MOs spectrum.^ 

Nevertheless, new types of molecular solids have been re- 
cently synthesized, for instance, aiming at the metalization 
of the single-component molecular solids family M(tmdt)2 
(M = Ni, Au),^ where tmdt stands for trimethylenetetrathia- 
fulvalenedithiolate. The molecular extension of the tmdt sys- 
tem is so large that a description based on a single molec- 
ular orbital (MO) is questionable. As a matter of fact, 
ab initio calculations have been performed^^^^ and an ef- 
fective three-band Hubbard model has been proposed and 
succeeded in describing electronic structures. In this re- 
spect, the quasi-one-dimensional molecular compound (TTM- 
TTP)l3/'^-^^ where TTM-TTP=2,5-bis(4,5-bis(methylthio)- 
l,3-dithiol-2-ylidene)-l,3,4,6-tetrathiapentalene looks like a 
promising candidate to investigate the theory limitation of a 
single-MO approximation. The presence of an I3 counter- 
anion is suggestive of formally organic [TTM-TTP] + cations 
in the crystal structure, i.e., unpaired electrons localized 
within the organic moieties. The magnetic susceptibility ^^'^^ 
and NMR measurements ^^'^^ revealed phase transitions at fi- 
nite temperature. Insulating and non-magnetic behaviors have 
been confirmed at low- temperature. Furthermore, based on 



Raman- scattering^ and x-ray measurements,^"^ it has been 
suggested that an asymmetric deformation of the [TTM-TTP] 
molecule occurs and charge disproportionation within the 
molecule is possible below a transition temperature. This 
novel charge-ordered (CO) state is different from the con- 
ventional CO state^ and is called intramolecular CO state. 
This nontrivial phenomena observed in this particular com- 
pound may not be described using a single-band description, 
i.e., the present molecular assembly could give rise to a phe- 
nomenon beyond a concept expected from a knowledge of 
a single molecule although the trigger may be hidden in the 
properties of multi molecular orbitals. Thus, there is a cru- 
cial need for theoretical description to rationalize the origin 
of such state. The purposes of our study are (i) to show that 
the energy of singly occupied molecular orbital (SOMO) and 
that of second-highest-occupied molecular orbital (HOMO- 
1) are quasidegenerate, (ii) to clarify the origin of this quasi- 
degeneracy using an effective three-fragment model for the 
[TTM-TTP] molecule, and (iii) to show explicitly that the in- 
tramolecular charge ordering actually occurs in the neighbor- 
ing two-molecular system. With this goal in mind, multirefer- 
ence wave function-based ab initio calculations are performed 
to rationalize the electronic distribution in the (TTM-TTP)l3 
material. 



II. MOs OF THE [TTM-TTP1+ ION 

The [TTM-TTP] organic molecule C14S12H12 is shown in 
Fig. [U The atomic coordinates are read from the (TTM- 
TTP)l3 298 K crystal structure. In order to assess the pack- 
ing influence in the electron trapping phenomenon, we did not 
perform any geometry optimization. Based on this structure, 
we performed correlated ab initio calculations. This type of 
approach is very insightful since important information is ac- 
cessible through a reading of the wave function. In particular, 
complete active space self-consistent field (CASSCF) calcula- 
tions have turned out to be very efficient to unravel intriguing 
electronic distribution in organic radical-based materials.—"— 
On top of the CASSCF wave functions, second-order pertur- 
bation theory calculations (CASPT2) were performed using 
an imaginary shift of 0.3 a.u. and an ionization potential- 
electronic affinity (IPEA) shift of 0.0 a.u. The IPEA shift 
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FIG. 1: [TTM-TTP] molecule (right-top) and energy levels of the 
isolated [TTM-TTP]^ ion from ROHF ab initio SCF calculations. 
The SOMO (u) and HOMO-l (g) are drawn where the red and blue 
colors denote the sign of orbital. 




FIG. 2: Effective three-fragment model for the u MO (a) and the g 
MO (b). The three fragments are named as L, R, and C, represent- 
ing the left, right, and center fragments, and their wave functions are 
^L, v^R, and cpc, respectively. The u MO is given by the superposi- 
tion of (/?L, — ^c, and cpR, while the g MO is made of cpi, and —(pR. 
The mutual interaction between the respective fragments are th and 
ts, representing the through-bond and through- space interactions, re- 
spectively. 



aims at correcting the energy differences calculations between 
states holding different open shells. Since we are dealing with 
spin states, it would be irrelevant to turn on this parameter. ^^'^^ 
This procedure allows one to incorporate the important dy- 
namical correlation effects to reach a high level of accuracy. 
All our ab initio calculations were performed using the MOL- 
CAS package^ ^ with all electron basis sets contractions for 
the elements S(7s6pld)/[4s3pld], C(5s5pld)/[3s2pld] and H 
(3s)/[ls]. We checked the validity of these particular contrac- 
tions by including diffuse and polarization functions which 
did not lead to any quantitative changes. 

At room- temperature, the [TTM-TTP] molecule exhibits an 
inversion center. Thus, the IVlOs can be classified as ger- 
ade (g) or ungerade (u) according to the symmetry point 
group. As mentioned before, the (TT]V[-TTP)l3 material is 
a charge-transfer salt, consisting of [TT]V[-TTP]+ and I3 
species. In a simple picture, the HOIMO of the [TT]V[-TTP]+ 
ion is half-filled, i.e., SOIMO. It has been usually recognized 
that one may concentrate on this SOIVEO, ignoring the rest 
of the spectrum on the assumption that its energy is well 
isolated from those of the other orbitals as compared to the 
bandwidth. This is one particular issue we wanted to exam- 
ine. Thus, semiempirical extended Hiickel calculations were 
first performed.^2'^^ In the following, the SOIMO and HOIMO- 
1 will be referred to as the u and g orbitals, respectively (see 
Fig. [T]). The energy separation between the g and u valence 
IVlOs (see Fig. [B is ~ 0.2 eV, while the bandwidth of the 
HOIMO is ?^ 1 eV.^^ Thus, the effect of the g orbital might not 
be negligible since it is likely to participate in the intramolecu- 
lar CO phenomenon expected in the (TT]V[-TTP)l3 system. In 
order to clarify the charge distribution, CAS SCF calculations 
were then carried out allowing the occupation of two IMOs 
by three electrons, i.e., CAS [3,2]. This method is known to 
provide very satisfactory charge distribution as soon as the 
active space is flexible enough. The g and u IMOs are then 
treated on the same footing and both symmetries states can 
be examined along these calculations. The CAS SCF energy 
difference between the g and u doublets is ^ 0.5 eV, confirm- 
ing the relative proximity of the frontier orbitals. The energy 
spectrum was finally calculated using a restricted open- shell 



SCF (ROHF) procedure. Along these ROHF calculations, 
three electrons are likely to occupy the frontier IMOs u and 
g. These IMOs are likely to be singly occupied or doubly oc- 
cupied in the CAS SCF calculations. In order to position the 
corresponding energies, we performed a ROHF calculation as- 
suming an average number of electrons in the g and u IMOs, 
i.e., 1.5 electrons. The calculated energy levels in the vicin- 
ity of the SOIMO are shown in Fig. [T] From this inspection, 
the SOIMO has ungerade character, whereas the HOlMO-l is 
gerade type in agreement with our extended Hiickel calcula- 
tions. The respective ROHF energies = — 8.27 eV and 
Eg = —8.69 eV, while the energy difference between the u and 
g IMOs is ~ 0.4 eV, a value which is consistent with our ex- 
tended Hiickel estimation. This combined semiempirical and 
ab initio information upon the constitutive unit [TT]V[-TTP]+ 
raises the relevance of a one-band approach to examine the 
electronic properties of (TT]V[-TTP)l3 crystal. 



III. EFFECTIVE THREE-FRAGMENT MODEL 

In this section, we analyze the chemical origin of the close- 
in-energy character of the g and u IMOs. Let us split the [TTIVL- 
TTP] molecule into 3 fragments — L, R, and C — as shown in 
Fig. [2I The L and R parts represent the left and right parts 
of the [TTIVL-TTP] molecule, whereas C corresponds to the 
ethylene-type bridging moiety. Based on this fragments pic- 
ture, the u MO displays "bonding" character while the g MO 
is "anti-bonding". Therefore, one may wonder why the bond- 
ing MO lies higher in energy than the anti-bonding one. Part 
of the answer can be found in the C group orbital u which is 
likely to mix in the L-R "bonding" orbital (see Fig. As 
for the interaction between these fragments, we consider two 
types of hopping integrals here, namely, and 4- The for- 
mer accounts for the through-bond interaction, while the latter 
represents the through- space interaction. The energy levels of 
the isolated L and R fragments are identical, and set to Eq. On 
the other hand, the energy level of the isolated C fragment is 
much lower than the L and R fragments, and we parametrize 
the energy difference as Ae (> 0). Using the local orbitals 



(/?L, ^R, and (/?c shown in Fig. O the effective Hamiltonian 
for the 3 -fragment model reads 



-ts So -th R 
-tb -^b eo- Ae I c 



(1) 



First, let us neglect the through-bond th interaction. The 
through- space interaction between the L and R parts results 
in bonding and antibonding MOs whose energy difference is 
2ts. This virtual situation is shown in the center of Fig. [3] As 
soon as the interaction is turned on, the energy of the bond- 
ing MO is shifted higher, while the anti-bonding MO remain 
unchanged due to symmetry constraints. Assuming that Ae is 
much larger than | tb | , the eigenvalues read 
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The respective wave functions are given by 
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where the quantity a is given hy a ^ (2t\^/ Ae). In this pic- 
ture, the L, R and C fragments wave functions are assumed 
to be orthogonal which is obviously not a limitation for the 
description. By comparing Eq. © with the u and g MOs 
shown in Fig. [TJ one can conclude distinctly that the wave 
functions Lp^ and Lp^ correspond to the u and g MOs, respec- 
tively. To scope out the Lp^' MO which is given in Eq. (l3cl) . 




FIG. 4: The L and R MOs obtained through a unitary transformation 
of Eq. 0]) by taking a = 0.25. 



we look into the set of ungerade MOs obtained from the SCF 
calculations. The in-phase combination uinner of Lpi,, (/^r, 
and (pc wave functions is shown in Fig [3] and its energy is 
^Uinner = —19.59 cV. Slncc di typical value for |tb| is of the 
order of several eV, the latter energy separation fully justifies 
the approximate expressions given previously. Thus, we can 
assign the (p^f wave function to the uinner MO identified in the 
SCF calculations. 

The ratio 2t\^/ Ae{'^ a) was estimated from the charge 
distribution in the [TTM-TTP]+ ion, i.e., so-called Mul- 
liken charges. Based upon the CAS [3,2] wave functions, 
CASPT2 calculations were performed and the resulting Mul- 
liken charges were estimated. Following the three-fragment 
model, we summed up the Mulliken charges within each frag- 
ment, e.g., pi, = Xl^^L P'i where pi is the charge of the re- 
spective atoms. In the ground state configuration where the u 
MO is singly occupied, the Mulliken charges on the respective 
fragments of the [TTM-TTP]+ ion were estimated as 



PL = PR = +0.634, PC 



-0.268. 



(4) 



The total charge is pl + Pr + Pc = +1- In order to ex- 
tract the contribution from the SOMO, we analyzed the charge 
difference between the cation [TTM-TTP]+ and the [TTM- 
TTP] molecule. The latter exhibits the following Mulliken 
charges on the respective fragments Pl = Pr = +0.149 and 
Pq = —0.297. Thus, the charge differences A = {pz—p%) 
(Z = L, R, C) were calculated as 



ApL ApR 0.485, Ape 0.029. 



(5) 
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FIG. 3: MOs diagram for the effective three-fragment model of the 
[TTM-TTP] molecule. As the through-space interaction U is turned 
on, the L and R local orbitals mix in and result in bonding and anti- 
bonding MOs (center). Then, the through-bond interaction th pushes 
the bonding MOs higher in energy (right). The identified MOs, u, g, 
and, Uinner, wMch are obtained from the SCF calculations, are also 
shown. The respective eigen-energies are given by = —8.27 eV, 
= -8.69 eV, and £u.._. = -19.59 eV. 



This means that the hole is mainly localized on the L and 
R fragments, whereas the charge on the C fragment remains 
almost unchanged. By combining the information upon the 
SOMO [Eq. ([3al) 1 and the numerical data of the charge distri- 
bution [Eq. Q], we can estimate the parameter a ^ (2tb / Ae) 
by using the relations 1/(2 + a^) = ApL or a^/(2 + a^) = 
Apc. Based on this analysis, we find a ^ 0.25 and finally 
the ratio (tb/Ae) is ^ 0.12. Once the parameter a is deter- 
mined, the fragment wave functions can be obtained explicitly 
through the inverse unitary transformation of Eq. ([3]). The re- 
sulting L and R MOs are shown in Fig. |4l and correspond to 
the and (fR local MOs. Since [{e^ + eg)/2 - ^Uinner] = 
Ae[l + 3(tb/A£)2] [see Eq. 0], the ROHF eigenvalues 
eu = -8.27 eV, e^ = -8.69 eV, and skinner = -19.59 eV 
lead to the effective energy of the C fragment Ae ^ 10.6 eV. 
Finally, the magnitudes of the through-bond/through- space in- 
teractions can be estimated, tb ^ 1-32 eV and t^ ~ —0.04 eV. 
From our evaluation, t^ is almost negligible, and the bonding 
and antibonding MOs of Fig. [5] are almost degenerated. The 
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through- space integral is much smaller than the through-bond 
one, |tb/ts| ~ 28. Since the energy difference between the u 
and g MOs is determined from Eq. ^ 

£u-%^-2ts + ^, (6) 

the origin of the quasidegeneracy of the u and g MOs can be 
clarified in the light of the extracted parameters. 

In conclusion, the u and g MOs ordering is completely 
determined by the relative energetics of this three-piece 
molecule. The through-space interaction ts is almost negligi- 
ble, while through-bond interaction determines the energy 
spectrum. 

IV. INTRAMOLECULAR CHARGE ORDERING 




FIG. 5: Valence MOs of the [TTM-TTP]2+ system obtained from 
CAS[6,4]SCF ab initio calculations upon the singlet state of g sym- 
metry. 



In this section, the interactions between neighboring [TTM- 
TTP]+ cations have been investigated using a dimer extracted 
from the (TTM-TTP)l3 crystal. We used the atomic coordi- 
nates of the crystal structure at 298 K. If the crystal- structure 
data at low-temperature phase were available, we may per- 
form more quantitative analysis by combining the geometrical 
optimization. However, it can be considered that the following 
results are not affected qualitatively. Since we were not only 
interested in charge distribution, CASPT2 calculations were 
also performed to specify the low-energy spectroscopy of the 
[TTM-TTP]^+ dimer. 

A. "MOs" 

First, CASSCF calculations were carried out on the [TTM- 
TTP]2^ dimer. These calculations were performed including 
six electrons in four MOs in the active space (CAS [6,4]) to 
account for the important static correlation effects. The inter- 
actions between the neighboring [TTM-TTP]+ ions labeled as 

1 and II give rise to the effective MOs shown in Fig. [5] Impor- 
tantly, the dimer exhibits an inversion center and the resulting 
MOs can be classified into gerade and ungerade. Neverthe- 
less, the inversion center of each individual subunit is lost, 
which might lead to electron localization within the [TTM- 
TTP]+ building blocks. As expected, the CAS[6,4]SCF fron- 
tier orbitals (Fig. [5]) consist in the in-phase and out-of-phase 
combinations of the (pi^^i, (pR^i, V^L,n, and (pR^ji MOs. The 
(/?L,i and (/9rj are the left (L) and right (R) localized MOs on 
the subunit I (see Fig. |4]). A similar definition holds for the 
subunit II. 

From the subsequent CASPT2 treatment, the ground state is 
a gerade singlet, and the active MOs exhibit occupation num- 
bers 2.0, 2.0, 0.9, and 1.1, respectively. Importantly, the oc- 
cupation numbers of the g2 and U2 MOs strongly deviate from 

2 and 0, a feature of the open-shell nature of the ground state 
singlet. The symmetry breaking within each [TTM-TTP]+ 
ion induced by the presence of a second partner leads to a 
mixing of the orbitals depicted in Fig. [3] This particular mech- 
anism is likely to result in the charge ordering we now wish to 
examine. 



B. Mulliken charge 

From the CASPT2 calculations performed upon the [TTM- 
TTP]2^ dimer, a similar Mulliken charge analysis was car- 
ried out to quantify the charge redistribution accompanying 
the dimer formation. The ground singlet state exhibits the 
following charge reorganization as compared to the isolated 
neutral [TTM-TTP] molecule 

ApL = 0.697, ApR = 0.272, Apc = 0.030. (7) 

As a major conclusion, the charge on the left (right) fragment 
is strongly enhanced (reduced), which is a clear indication of 
the intramolecular charge ordering. Again, the charge on the 
center fragment C is almost unchanged. The present eval- 
uation based on the two-molecule dimer overestimates the 
charge difference between the L and R fragments as com- 
pared to its value in the crystal, and also the validity of Mul- 
liken charges remains questionable since the fluctuation ef- 
fects were not included. Nevertheless, this qualitative analysis 
supports the intramolecular CO state. 

C. Low-energy spectroscopy 

Our Mulliken charge analysis is suggestive of an in- 
tramolecular charge ordering mechanism in the (TTM-TTP)l3 
material. Thus, one may expect a charge localization on the 
left or right parts of the TTM-TTP building blocks. Such sce- 
nario is likely to give rise to magnetic interactions involving 
either the inner parts (i.e., (/prj and (/:?L,n fragments MOs) or 
the outer parts (i.e., and (fR^ji). 

Thus, starting from the CAS[6,4]SCF calculations, the low- 
energy spectroscopy of the [TTM-TTP] 2^ species was in- 
spected. The multi-reference CASSCF wave functions were 
expanded using a local orbitals basis set {(ph,i, ^r,i, ^L,n, 
^R,n} following the transformation 

^gi ^ -^(-^R,i + ^L,n), (8a) 

^ui ~ ^(-(/9R,i - (/;L,n), (8b) 
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TABLE I: Wave functions decompositions (weights) of the [TTM-TTP]2^ low-energy states. The wave functions are expressed in terms of 
the two hole determinants. CAS[6,4]PT2 energies with respect to the ground state singlet g are given in eV. 

|v^L,i^^R,n) |v^R,i^^L,n) I ^/?L,i ^^L,n) |v^L,i^^L,i) | V^R,i ^^R,i) Energy 

(%) (%) (%) (%) (%) (%) (eV) 

Singlet, gerade 84 1 3 9 3 Eq = 0.00 

n 27 44 13 5 E2 = 0.27 

\(pL,i(pR,Ji) \(pR,i(pL,Ji) \(pL,i(pR,i) \(pL,i(pL,Ji) Energy 

(%) (%) (%) (%) (eV) 

Triplet, ungerade 62 5 33 ^1 = 0.09 

28 47 2 23 E3 = 0.47 

6 8 70 16 ^4 = 0.62 



^g2 ^ --^(+(/?L,i - ^R,n), (8c) 
^u2 ~ --^(-(/?L,i - ^R,n). (8d) 

Since we are dealing with a six-electron/four-MO system, 
we performed our analysis in the two-hole/four-MO picture. 
Schematic representations of the different hole configurations 
for the singlet and triplet states based on these local orbitals 
are shown in Fig.[6l Such transformation affords a reading of 
the different wave functions and the extractions of the relevant 
information in a valence-bond type analysis. Since | is rel- 
atively small, the left and right moieties within each unit do 
not significantly overlap. Thus, the previous transformation is 
almost unitary. 

The lowest lying states of the [TTM-TTP]2^ species are 
given in Table U The ground state is singlet g and consis- 
tent with a picture that minimizes the electrostatic energy be- 
tween the two [TTM-TTP]+ ions. The comparison between 
the ground state singlet g and first triplet u is very instruc- 
tive. In both cases, the wave function is largely dominated by 
the electronic configurations (84 % and 62 % in Table U for 
the singlet and triplet states, respectively) involving the inner 
parts (see Fig. [6]). The energy difference affords an evaluation 
of the intermolecular magnetic exchange interaction : 

J = Ei-Eo^ 0.09 eV (9) 

A positive value reflects an antiferromagnetic interaction, re- 




I'^mV^E.n) I^Rj^L,n> IV^LJ^R,l> I^L,i^L,n> I^LJ^Lj) I'^ej^'rj) 



lv^L.i^R,n> l¥^Rj¥'L,n> Iv^mV^Rj) lv^Lj¥'L,n> 

FIG. 6: Schematic view of the different configurations using a hole 
picture for the singlet (top) and triplet (bottom), respectively. The 
subscripts I and 11 stand for the molecules I and n as depicted in Fig. 

m 



suiting from the intramolecular charge ordering. Such a value 
suggests that the spin gap excitation energy should be ^ 900 
K. The first-excited singlet corresponds to the holes localiza- 
tion within one [TTM-TTP]+ ion. This is a reflection of an 
intermolecular charge ordering which lies much higher and 
might not be relevant to describe the low-energy properties of 
the (TTM-TTP)l3 material. The energy difference (0.35 eV) 
between this second excited state (singlet g) and the fourth 
one (triplet u) supports the strong antiferromagnetic interac- 
tion within each subunit. This is a mechanism that involves 
again intermolecular electron transfer. The low-energy spec- 
troscopy in the (TTM-TTP)l3 material should be mostly gov- 
erned by the through-space interdimer magnetic interaction J, 
which results from the intramolecular charge ordering. 

V. SUMMARY AND DISCUSSIONS 

In the present paper, we examined the low-energy prop- 
erties of the (TTM-TTP)l3 material. Complementary semi- 
empirical and wave function-based ab initio calculations were 
performed upon the elementary unit and a dimer to investi- 
gate the underlying electronic distribution. We showed that 
the SOMO and HOMO-1 are close in energy and that a one- 
band picture should be ruled out. A chemical understanding 
arises from the inspection of the building block and effective 
parameters extracted from a three-fragment model allowed us 
to rationalize the quasidegeneracy of the [TTM-TTP] fron- 
tier MOs. From the calculation of a system consisting of 
two neighboring molecules, the inversion symmetry within 
each molecule is lost. This result is in agreement with ex- 
perimental data. Indeed, the intramolecular charge ordering 
which has been experimentally suggested actually occurs in 
the two-molecule cluster system. The mixing of the SOMO 
and HOMO-1 results in an electron trapping which results in 
a 0.697 / 0.272 on the L / R moieties of the [TTM-TTP] + 
ion. Finally, our ab initio calculations suggest that the low- 
energy of the (TTM-TTP)l3 material is controlled by a single 
exchange interaction 0.09 eV, resulting from the charge order- 
ing. 

In order to complement our Mulliken charge analysis, we 
also performed preliminary Raman and infrared (IR) calcu- 
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lations upon the [TTM-TTP]+ unit. It has been observed 
upon cooHng that the Raman 1490 cm~^ band spHts into two 
peaks centered at 1487 cm~^ and 1499 cm~^, where this 
phenomenon was attributed to the differentiation between the 
C=C yhdene bonds, featuring a symmetry breaking. We 
used the GAUSSIAN03 package^^ on isolated one- and two- 
molecule systems, and full geometry optimizations were car- 
ried out on both systems. The information extracted from 
the one-molecule system should be compared with the high- 
temperature regime. In contrast, the two-molecule system 
is expected to give access to the low-temperature Raman 
spectrum characteristics. The calculated spectrum of two- 
molecule system displays two vibrational frequencies, 1522 
cm~^ and 1534 cm~^, while that of the 1 -molecule system 
shows a single band at 1535 cm~^. Despite a general blue- 
shift, this result is in agreement with experimental data, and 
can be attributed to an intramolecular CO phase accompany- 
ing the descent in symmetry within the TTM-TTP units. In 
this calculation, the [TTM-TTP] 2^ dimer was extracted from 
the available crystallographic data. A direct comparison with 
experimental findings would rely on Raman and IR calcula- 
tions using the low-temperature crystal structure. Unfortu- 



nately, the lack of such x-ray data disposes of this strategy. 

Throughout this work, we have focused only on the [TTM- 
TTP] ion and have neglected the effect of the counteranion 
of I3 . As a matter of fact, it has been pointed out, from 
X-ray measurement,^^ that the large displacement of the I3 
species occurs in the low-temperature phase. This might yield 
changes in the electrostatic potential, and possibly the elec- 
tronic distribution in the [TTM-TTP] ion would qualitatively 
be modified. This point is to be studied in future work. 
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